The database of experimentally supported targets: a functional update of TarBase

TarBase5.0 is a database which houses a manually curated collection of experimentally supported microRNA (miRNA) targets in several animal species of central scientific interest, plants and viruses. MiRNAs are small non-coding RNA molecules that exhibit an inhibitory effect on gene expression, interfering with the stability and translational efficiency of the targeted mature messenger RNAs. Even though several computational programs exist to predict miRNA targets, there is a need for a comprehensive collection and description of miRNA targets with experimental support. Here we introduce a substantially extended version of this resource. The current version includes more than 1300 experimentally supported targets. Each target site is described by the miRNA that binds it, the gene in which it occurs, the nature of the experiments that were conducted to test it, the sufficiency of the site to induce translational repression and/or cleavage, and the paper from which all these data were extracted. Additionally, the database is functionally linked to several other relevant and useful databases such as Ensembl, Hugo, UCSC and SwissProt. The TarBase5.0 database can be queried or downloaded from http://microrna.gr/tarbase

Functional Transcriptomics in Diverse Intestinal Epithelial Cell Types Reveals Robust MicroRNA Sensitivity in Intestinal Stem Cells to Microbial Status

Gut microbiota play an important role in regulating the development of the host immune system, metabolic rate, and at times, disease pathogenesis. The factors and mechanisms that mediate interactions between microbiota and the intestinal epithelium are not fully understood. We provide novel evidence that microbiota may control intestinal epithelial stem cell (IESC) proliferation in part through microRNAs (miRNAs). We demonstrate that miRNA profiles differ dramatically across functionally distinct cell types of the mouse jejunal intestinal epithelium and that miRNAs respond to microbiota in a highly cell type-specific manner. Importantly, we also show that miRNAs in IESCs are more prominently regulated by microbiota compared with miRNAs in any other intestinal epithelial cell subtype. We identify miR-375 as one miRNA that is significantly suppressed by the presence of microbiota in IESCs. Using a novel method to knockdown gene and miRNA expression ex vivo enteroids, we demonstrate that we can knock down gene expression in Lgr5+ IESCs. Furthermore, when we knock down miR-375 in IESCs, we observe significantly increased proliferative capacity. Understanding the mechanisms by which microbiota regulate miRNA expression in IESCs and other intestinal epithelial cell subtypes will elucidate a critical molecular network that controls intestinal homeostasis and, given the heightened interest in miRNA-based therapies, may offer novel therapeutic strategies in the treatment of gastrointestinal diseases associated with altered IESC function

miR-30 Family Controls Proliferation and Differentiation of Intestinal Epithelial Cell Models by Directing a Broad Gene Expression Program That Includes SOX9 and the Ubiquitin Ligase Pathway

Proliferation and differentiation of intestinal epithelial cells (IECs) occur in part through precise regulation of key transcription factors, such as SOX9. MicroRNAs (miRNAs) have emerged as prominent fine-tuners of transcription factor expression and activity. We hypothesized that miRNAs, in part through the regulation of SOX9, may mediate IEC homeostasis. Bioinformatic analyses of the SOX9 3′-UTR revealed highly conserved target sites for nine different miRNAs. Of these, only the miR-30 family members were both robustly and variably expressed across functionally distinct cell types of the murine jejunal epithelium. Inhibition of miR-30 using complementary locked nucleic acids (LNA30bcd) in both human IECs and human colorectal adenocarcinoma-derived Caco-2 cells resulted in significant up-regulation of SOX9 mRNA but, interestingly, significant down-regulation of SOX9 protein. To gain mechanistic insight into this non-intuitive finding, we performed RNA sequencing on LNA30bcd-treated human IECs and found 2440 significantly increased genes and 2651 significantly decreased genes across three time points. The up-regulated genes are highly enriched for both predicted miR-30 targets, as well as genes in the ubiquitin-proteasome pathway. Chemical suppression of the proteasome rescued the effect of LNA30bcd on SOX9 protein levels, indicating that the regulation of SOX9 protein by miR-30 is largely indirect through the proteasome pathway. Inhibition of the miR-30 family led to significantly reduced IEC proliferation and a dramatic increase in markers of enterocyte differentiation. This in-depth analysis of a complex miRNA regulatory program in intestinal epithelial cell models provides novel evidence that the miR-30 family likely plays an important role in IEC homeostasis

Genome-Wide Analysis of Natural Selection on Human Cis-Elements

Background: It has been speculated that the polymorphisms in the non-coding portion of the human genome underlie much of the phenotypic variability among humans and between humans and other primates. If so, these genomic regions may be undergoing rapid evolutionary change, due in part to natural selection. However, the non-coding region is a heterogeneous mix of functional and non-functional regions. Furthermore, the functional regions are comprised of a variety of different types of elements, each under potentially different selection regimes. Findings and Conclusions: Using the HapMap and Perlegen polymorphism data that map to a stringent set of putative binding sites in human proximal promoters, we apply the Derived Allele Frequency distribution test of neutrality to provide evidence that many human-specific and primate-specific binding sites are likely evolving under positive selection. We also discuss inherent limitations of publicly available human SNP datasets that complicate the inference of selection pressures. Finally, we show that the genes whose proximal binding sites contain high frequency derived alleles are enriched for positive regulation of protein metabolism and developmental processes. Thus our genome-scale investigation provide

Transcriptional targets of the schizophrenia risk gene MIR137

Genome-wide association studies (GWAS) have strongly implicated MIR137 (the gene encoding the microRNA miR-137) in schizophrenia. A parsimonious hypothesis is that a pathway regulated by miR-137 is important in the etiology of schizophrenia. Full evaluation of this hypothesis requires more definitive knowledge about biological targets of miR-137, which is currently lacking. Our goals were to expand knowledge of the biology of miR-137 by identifying its empirical targets, and to test whether the resulting lists of direct and indirect targets were enriched for genes and pathways involved in risk for schizophrenia. We overexpressed miR-137 in a human neural stem cell line and analyzed gene expression changes at 24 and 48 h using RNA sequencing. Following miR-137 overexpression, 202 and 428 genes were differentially expressed after 24 and 48 h. Genes differentially expressed at 24 h were enriched for transcription factors and cell cycle genes, and differential expression at 48 h affected a wider variety of pathways. Pathways implicated in schizophrenia were upregulated in the 48 h findings (major histocompatibility complex, synapses, FMRP interacting RNAs and calcium channels). Critically, differentially expressed genes at 48 h were enriched for smaller association P-values in the largest published schizophrenia GWAS. This work provides empirical support for a role of miR-137 in the etiology of schizophrenia

Differential impact of glucose administered intravenously and orally on circulating mir-375 levels in human subjects

Background: To date, numerous nucleic acid species have been detected in the systemic circulation including microRNAs (miRNAs); however their functional role in this compartment remains unclear. Objective: The aim of this study was to determine whether systemic levels of miRNAs abundant in blood, including the neuroendocrine tissue-enriched miR-375, are altered in response to a glucose challenge. Design: Twelve healthy males were recruited for an acute cross-over study which consisted of two tests each following an eight-hour fasting period. An oral glucose tolerance test (OGTT) was performed and blood samples were collected over a 3-hour period. Following a period of at least one week, the same participants were administered an isoglycemic intravenous glucose infusion (IIGI) with the same blood collection protocol. Results: The glucose response curve following the IIGI mimicked that obtained after the OGTT, but as expected systemic insulin levels were lower during the IIGI compared to the OGTT (P<0.05). MiR-375 levels in circulation were increased only in response to an OGTT and not during an IIGI. In addition, the response to the OGTT also coincided with the transient increase of circulating glucagon-like peptide-1 (GLP-1), glucagon-like peptide-2 (GLP-2), and glucose-dependent insulinotropic polypeptide (GIP). Conclusions: The present findings show levels of miR-375 increase following administration of an OGTT and in light of its enrichment in cells of the gut, suggest that the gastrointestinal tract may play a significant role to the abundance and function of this microRNA in the blood

Comprehensive analysis of The Cancer Genome Atlas reveals a unique gene and non-coding RNA signature of fibrolamellar carcinoma

Fibrolamellar carcinoma (FLC) is a unique liver cancer primarily affecting young adults and characterized by a fusion event between DNAJB1 and PRKACA. By analyzing RNA-sequencing data from The Cancer Genome Atlas (TCGA) for >9,100 tumors across ~30 cancer types, we show that the DNAJB1-PRKACA fusion is specific to FLCs. We demonstrate that FLC tumors (n = 6) exhibit distinct messenger RNA (mRNA) and long intergenic non-coding RNA (lincRNA) profiles compared to hepatocellular carcinoma (n = 263) and cholangiocarcinoma (n = 36), the two most common liver cancers. We also identify a set of mRNAs (n = 16) and lincRNAs (n = 4), including LINC00473, that distinguish FLC from ~25 other liver and non-liver cancer types. We confirm this unique FLC signature by analysis of two independent FLC cohorts (n = 20 and 34). Lastly, we validate the overexpression of one specific gene in the FLC signature, carbonic anhydrase XII (CA12), at the protein level by western blot and immunohistochemistry. Both the mRNA and lincRNA signatures support a major role for protein kinase A (PKA) signaling in shaping the FLC gene expression landscape, and present novel candidate FLC oncogenes that merit further investigation

Accurate microRNA target prediction correlates with protein repression levels

MicroRNAs are small endogenously expressed non-coding RNA molecules that regulate target gene expression through translation repression or messenger RNA degradation. MicroRNA regulation is performed through pairing of the microRNA to sites in the messenger RNA of protein coding genes. Since experimental identification of miRNA target genes poses difficulties, computational microRNA target prediction is one of the key means in deciphering the role of microRNAs in development and diseas

A Systematic Mapping Approach of 16q12.2/FTO and BMI in More Than 20,000 African Americans Narrows in on the Underlying Functional Variation: Results from the Population Architecture using Genomics and Epidemiology (PAGE) Study

Genetic variants in intron 1 of the fat mass- and obesity-associated (FTO) gene have been consistently associated with body mass index (BMI) in Europeans. However, follow-up studies in African Americans (AA) have shown no support for some of the most consistently BMI-associated FTO index single nucleotide polymorphisms (SNPs). This is most likely explained by different race-specific linkage disequilibrium (LD) patterns and lower correlation overall in AA, which provides the opportunity to fine-map this region and narrow in on the functional variant. To comprehensively explore the 16q12.2/FTO locus and to search for second independent signals in the broader region, we fine-mapped a 646-kb region, encompassing the large FTO gene and the flanking gene RPGRIP1L by investigating a total of 3,756 variants (1,529 genotyped and 2,227 imputed variants) in 20,488 AAs across five studies. We observed associations between BMI and variants in the known FTO intron 1 locus: the SNP with the most significant p-value, rs56137030 (8.3×10-6) had not been highlighted in previous studies. While rs56137030was correlated at r2>0.5 with 103 SNPs in Europeans (including the GWAS index SNPs), this number was reduced to 28 SNPs in AA. Among rs56137030 and the 28 correlated SNPs, six were located within candidate intronic regulatory elements, including rs1421085, for which we predicted allele-specific binding affinity for the transcription factor CUX1, which has recently been implicated in the regulation of FTO. We did not find strong evidence for a second independent signal in the broader region. In summary, this large fine-mapping study in AA has substantially reduced the number of common alleles that are likely to be functional candidates of the known FTO locus. Importantly our study demonstrated that comprehensive fine-mapping in AA provides a powerful approach to narrow in on the functional candidate(s) underlying the initial GWAS findings in European populations

DIANA-microT web server: elucidating microRNA functions through target prediction

Computational microRNA (miRNA) target prediction is one of the key means for deciphering the role of miRNAs in development and disease. Here, we present the DIANA-microT web server as the user interface to the DIANA-microT 3.0 miRNA target prediction algorithm. The web server provides extensive information for predicted miRNA:target gene interactions with a user-friendly interface, providing extensive connectivity to online biological resources. Target gene and miRNA functions may be elucidated through automated bibliographic searches and functional information is accessible through Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. The web server offers links to nomenclature, sequence and protein databases, and users are facilitated by being able to search for targeted genes using different nomenclatures or functional features, such as the genes possible involvement in biological pathways. The target prediction algorithm supports parameters calculated individually for each miRNA:target gene interaction and provides a signal-to-noise ratio and a precision score that helps in the evaluation of the significance of the predicted results. Using a set of miRNA targets recently identified through the pSILAC method, the performance of several computational target prediction programs was assessed. DIANA-microT 3.0 achieved there with 66% the highest ratio of correctly predicted targets over all predicted targets. The DIANA-microT web server is freely available at www.microrna.gr/microT